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16. A6»itoci , . . i. j 

iSxperiinental flow InveBtlgationB on smooth airfoils were done 
using numerical solutions for transsonic airfoil streaming 
with shockless supersonic range. The experimental flow 
reproduced essential sections of the theoretically computed 
frictionless solution*. Agreement is better in the expansion 
part of the flow than in the compression part. The flow 
was nearly stationary in the entire Velocity range 
investigated. 
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The Studied Airfoil 


The starting point for our experiinental flow investigations 
on smooth airfoils was the numeric solutions for transsonic air- 
foil streaming with shockless supersonic range computed by G.Y. 
Nieuwland [1] and J.W. Boers toel [2] using the method of T.M. 
Cherry and M.J, Lighthill. Through these tests we shall determine 
whether this shockless compression is a physical possibility. 

The original symmetric airfoil provided for the experiments 
O.ilOO - 0.7500 - 1.3750 (fig. 1) was replaced after consultation 
with J.V). Boers toel of the NLR (National Lucht-en Ruimtevaart 
Laboratorium, Amsterdam) » by the profile 0.10 - 0.676 - 1.6 (fig. 
2), because the latter is supposed to have less tendency for 
parting of the laminar boundary layer at the airfoil nose and 
should give greater supersonic speeds at smaller Mach number of 
incident flow. This airfoil is still symmetric, it is 16.3% 
thick and is designed for an incident Mach number Ma = 0.745 with 
an incident-flow angle it = 0°. The greatest Mach number in the 
supersonic range — which is assumed to be far in front of the 
thickest point of the airfoil — amounts to Ma = 1.464. The profile 
coordinates, whose computation is very complex, were providad to 
us by the NLR. 

By using these coordinates, a profile template of ca. 635 mm 
length of chord was produced whose deviation from the theoretical 
value is less than 0.05 mm. This accuracy of application was 
easily achieved with the coordinatographs available to us. Using 
this template, a 5-fold smaller test model was produced by means 
of a duplicating milling machine. For a tight fit of spindles 
and joints in the machine, an accuracy of 0.01 mm can be attained 
in this manner. 

2 . Test Equipment 

The profile (fig. 3) is anchored with two each laterally 
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attached ire^*aining pins in appropriate drilled holes of the 35 mm 
thick canal window (fig. 4). This type of attachment prevents 
additional disturbances of the flow due to holding devices pro- 
truding into the canal. In addition, with this set-up, the 
set angle of the profile can be easily adjusted from the outside 
and measured. The depth of the profile is 100 mm (corresponding 
to canal depth), the length 126.8 mm and the max. thickness 20.68 
mm. The canal height is limited to 300 mm (fig. 5) by two canal 
walls slit in the flow direction; these walls can also be deformed, 
depending on their location (fig. 6). Thus, the flow generated in 
the canal can be adapted to the free profile streaming [3]. The 
opened canal rotated about a horizontal axis by 90° clearly shows 
Che half of the comb-like setting-dif fusor , with which the through- 
put is regulated. Figure 8 shows a complete picture of the test 
arrangement. The pipeline entering the picture on the upper right 
connects the test canal to the vacuum vessel via a fast-acting 
valve. Under the canal is a Mach-Zehnder Interferometer with 
attached high-frequency flash-lamp, to the left of this, the per- 
tinent power supply unit. Behind the Mach-Zehnder Interferometer 
we see a part of the differential interferometer. Both optic 
instruments are movable on rails and can be used optionally to 
view the streaming. 

The canal (wind tunnel) is operated intermittently with 
atmospheric air which flows into an evacuated vessel of 134 m 
after passing through the measurement lane. Feed to the canal 
with dry air was not possible for the tests described here since 
the pipeline from the existing air drier Co the canal is not 
properly sized for these experiments and thus an excessive pressure 
drop would result. At present, a gasometer of ca, 40 m volume 
is under development which will allow any adjustment of air humid- 
ity to permit a study of the influence of this parameter. 

For the investigation of flow with the Mach-Zehnder Interfero- 
meter, the instrument was adjusted for infinite strip width under 
a calm flow medium. In this case the interference strips correspond 
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to lines of constant density, if the deviations of the flow from 
the planar case can be neglected (e.g. wall boundary layer at the 
canal side-walls, secondary streaming). An estimation of the 
boundary layer thickness at the canal walls in the measured range 
(ignoring the Influence of the airfoil) gives a maximum value of 
0.02 mm. 

To register the interference pictures, a mirror-reflex 
camera and a high-speed camera (picture frequency up to 10 kHz) 
are used. A stroboscope is used as light source. The duration 
of the individual flashes is ca. 1 us, the converted electrical 
energy ca. 1 Joule. 

3. Study Results 

In general, we can Sky that the experimental flow reproduces 
essential sections of the theoretically computed frictionless 
solution, In the expansion part of the flow, the agreement is 
clearly better than in the compression part, since there the boun- 
dary layer increases quickly and its compression action changes 
the velocity field, The boundary layer is seen in the interfero- 
grams by the sudden bending vf Interference lines near the airfoil 
(e.g. fig. 9). It should be mentioned that the flow — except for 
the turbulent wake and the sound disruptions generated there — is 
nearly stationary in the entire velocity range investigated by us. 
One important prerequisite for this behavior lies in the large 
stability of the turbulent boundary layer in the range 0.2 < ^ 

against flow becoming non-laminar, even in the presence of compr*ess- 
ion shocks (e.g. fig. 13). 

*- 

The typioal run of the experimentally found flow state is 
presented in the interference pictures 9 - 13. The subsequent, 
dashed line gives the position of the theoretically calculated 
sound barrier for an incident flow Mach number M^q- 0.745. 

In our tests (rest pressure p^ = 745 mm Qs, rest temperature 20 *^0, 
depth of canal b » 100 mm, wavelength of used light 567 nm), 
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the sound barrier corresponds to interference strip with the 
order m * 17.25, if the order of the interference in the stagna- 
tion point of the airfoil is specified as m » 0 (see table). The 
Miich number of incident flow is determined from the interference 
pictures with an accuracy of 57e. In our previous investigations 
we found that in all cases in which a supersonic region appeared, 
a slanting compression shock goes out from the profile nose, which 
is reflected repeatedly between the sonic barrier and profile 
(under circumstances as an attenuation fan). The origin of the 
shock lies just behind the point of maximum curvc^ture of the profile 
nose (fig. 14). Since in the interference pictures the suspected 
non-laminar flow as source for the occurrence of the compression 
shock cannot definitely be found, the flow Was also examined by 
using the tinting method. For this purpose j the profile was 
subjected to a suspension of titanium dioxide in oil as spray, and 
Chen exposed ca* 5s to the flow. The tests showed that behind the 
point of maximum curvature, there is formation of a dead-water 
region of ca. 5 mm length. In the tinting pictures (fig. 19, 20) 
we see the position of the non-laminar region and that there the 
sprayed~on suspension remains on the profile, whereas both in front 
and behind it, the suspension has been blown off almost completely 
by Che flow. In addition, the pictures also show that the dead- 
water region does not extend up to the canal side walls, but ends 
ca. 10 mm in front of this. Apparently, the conversion of the 
laminal profile boundary layer occurs here, due to the effects of 
the turbulent wall boundary layer; it is earlier in this case so 
that the separation flow of the wall boundary layer does not 
occur at these places. 

As long as the observed supersonic range is larger or about 
equal to the calculated one, it is concluded by one or more weak, 
perpendicular, almost stationary compression shocks, whose origin 
we view as the flow field changed by the compression shock emanating 
from the profile nose. The disturbances emanating from the turbu- 
lent wake of the airfoil do not indicate any larger influence on 
the position and intensity of shocks (fig. 12, 13). If the super- 
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sonic range Is smaller than the calculated one, then it breaks down 
into a sequence of weak shocks and expansion fans. 

Another type of disturbance of the flow field is caused by 
the liberated heat on condensation of the water vapor contained 
in the atmospheric air. Figures 15-18 show the condensed water 
vapor in dark-field illumination. The relative humidity of the air 
was 50% as rest. The added heat can also be detected in the inter- 
ference strips. Figures 15, 16 correspond to flow state of fig. 11, 
figures 17, 18, correspond to that of fig. 13, 

Until now, no use had been made of the possibility mentioned 
above for adapting the canal contours (fig. 5, 6) to the flow lines 
expected in the free flow, since from other experiments it Was 
supposed that the longitudinal slit in the canal walls would 
cause the majority of the compensation flow. A sample change had 
no effect on the shape of the compression shocks at the profile 
front edge, whereas the other flow field was changed only to a 
minor degree. The position-dependent setting of the contours 
thus only comes into consideration for fine-tuning the flow outside 
the near-field of the profile. 

The pressure distributions calculated from the interferograms 

(fig. 9-15) along the profile are presented in figures 9a, b - 

13a, b. Instead of the pressure coefficient c in fig. 2, we 

Po * P P 

plotted the run of -5^ (p^ - stagnation pressure.^ p = static 

pressure) against the dimensionless length of chord since in 
our test set-up, the inflow velocity could not be determined very 
accurately from the interf erograros . However, the coordinate 
system has an affinity with that of fig. 2. As long as the 
bomidary layer is thin enough to be neglected, the pressure was 
determined from the order of interference at the particular pro- 
file location. As soon as the run of the interference strips 
indicates the shape of a boundary layer, the calculation of pressure 
takes place with the value of the interference order outside the 
boundary layer. 
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If we compare the obtained curves with those calculated from 
the theory, then we see a good qualitative agreement with regard 
to the profile, But on detailed consideration it turns out that 
due to the compression shock emanating from the "corner’* of the 
profile nose, attainment of the expected underpressure peak is 
prevented and thus the structure of the computed flow field is 
disrupted. For y > 0.2, the theoretical pressure distribution 
is falsified by the quickly increasing compression effect of the 
boundary layer. Even now we can say with certainly that for the 
region 0.2 4 ^ ^ 1 > the frictionless, calculated profile will 
have to be corrected by modification of the contour or suction- 
remova*;. in order to attain the calculated pressure distribution. 


6 


REFERENCES 


[Ij G.Y, Nt«dulaaclt Tiransontc potential flow around a faintly of 
quasi-elliptical aerofoil sections. NLR TR T 172) 1967, 

[2] J.W. Boerstoel: A survey of symmetrical transsonic potential 
flows around quasi-elliptical aerofoil sections. NLR TR T 
136, 1967. 

[3] Bernhard H. Goethert: Transonic Wind Tunnel Test Ing. 
AGARDograph 49, 1961, 


7 
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Kelation between interference etrip arrangement. Density, pressure, 
temperature and Mach-number. (Prerequisites; Isotropic flow, rest 
pressure, p ■ 745 mm QS, rest temp. T « 20 °C, light wavelength \m 
567 nm, depth of tunnel b ■ 100 mni) , ^ 


1,0 1,0 1,0 

0,989 0,905 0,99 

0,979 0.971 0.99 

0,9^8 0,956 0,98 

0,958 0.941 0.98 

0,947 0,925 0.97 

0,956 0.91 0,97 

0,926 0.896 0,96 

0,915 0,8820,96 

0,905 0,07 0,96 

0,894 0.8550,95 

0,885 0.64 0.95 

0,673 0.826 0,94 

0,862 0.811 0.94 

0.852 0,797 o.Sl 

0,841 0,783 0,93 

0.8304 0.771 0.9s 


0.756 

0.746 

0.735 

0.724 

0.714 

0.703 

0.693 

o:682 

0.671 

0.661 



Na 

1.0 

o.o5 

0.996 

0.15 

0.9916 

o.2o8 

0.9875 

0.25 

0.983 

0,295 

0.977 

0.333 

0.973 

0.37 

0,969 

0,390 

0.964 

0.424 

0.9605 

0.452 

0,956 

0,48 

0,951 

o.5o6 

0*946 

0.53 

0.9418 

0.555 

o.Sl? ' 

0.579 

0.932 

0.6 


0,62 


0,8092 0,742 0,918 0,66 

0,7986 0.73 0.9139 0.686 27.0 

0.788 0.715 0.908 0,71 27.5 

0.777 o.7o2 0.903 0.727 28.o 

0.7668 0.6BB 0.8975 o.75 


0.675 0.8925 o.77 
0,661 0,8875 o.79 
0,648 O.0625 0.81 
0.635 0.877 o.83 


0.693 0,597 0.862 
o:682 0,5830.857 
0.671 0.571 0,851 
0.661 0.56 0.846 
o.65o2 0.547 0.84 
0.6396 0.532 0.835 


17.5 0*6289 
18,0 0,6183 


18.5 

19.0 


23.5 


25.0 

25.5 


0.8198 0.755 0.9225 0.645 26. o 
0.8092 0,742 0,918 0.668 26.5 


0.623 0.0725 0.849 31. 0 
0.61 0.867 0,87 31.5 


0*6289 0.521 o.620 

0,6183 0,51 0,823 

0.6077 0.497 O.017 
0.5971 0,485 0.812 


0.89 32,0 0,3215 
0.91 32.5 0.3I09 
0,93 33.0 o,3oo3 
0.948 33.5 0.2897 
0.97 34.0 0.2791 
0.991 34.5 0.2685 


2o,o 0.5759 0.461 
2o.5 0.5653 o,45 


0.5653 0.45 0,794 

21.0 0.5547 0,438 0.788 

21.5 0.5441 0.426 0.782 

22.0 0.5335 0,416 0.777 
22*510.5229 o.4o3 0.769 

0.5123 0.39 0.763 

0.50I7 0.38 0.756 

0.4911 0.369 0.75 
0.4805 0.358 0.744 
0.4699 0.346 0.737 
0.4593 0.335 o,73 


0.335 0.73 
0,325 0.723 
0.314 0.717 
o,3o4 0.71 


0,4487 0,325 0.723 

0.4381 0.314 0.717 

0,4275 o,3o4 0,71 
0.4169 0.293 o,7o25 
0.4063 0.282 0.694 

28.5 0.3957 0.272 0.687 

29.0 0.3851 0.262 0.68 

29.5 0.3745 0.253 0.673 

30.0 0,3639 0.243 0,665 

30.5 0.3533 o.232 0.657 
0.3427 0.222 0,648 
0.3321 0.213 0.641 
0,3215 o.2o3 0,632 
0.3I09 0.194 0.624 
0.3003 0,185 0.615 
0,2897 0,175 o,6o5 
0.2791 0.167 0.597 

0,2685 0.15a 0,586 
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Figure 5; Tunnel Wall (Slit In Flow Direction) 



Figure 6: Wind-tunnel Wall (Variable Wall Shape) | 
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Fig. 9: Interferogram of the Flow for Mqq - 0.690 

■ calculated sonic barrier for Mqo ■ 0.745 



Fig. 10: » 0.73 
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Pressure Distribution on the Profile lower Side 
1 length of chord of the profile 
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Fig. 11a: Pressure Distribution on the Profile top Side 
1 length of chord of the profile 
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Fig. 11b; Pressure Distribution on the Profile lower Side 
1 length of chord of the profile 
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Fig. 13a: Pressure Distribution on the Profile top Side 
1 length of chord of the profile 



Fig. 13b: Pressure Distribution on the Profile lower Side 
1 length of chord of the profile 
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Figure 13: Water Vapor Condensation, M >0.78 
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Figure 16: Water Vapor Condensation, Mqq* 0.78 
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Detection of Boundary Layer separation by the tinting method 
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Figure 20: Lo*;er Side of the Profile 
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